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Abstract 



The excess of low invariant mass dilepton yield in heavy ion colli- 
sions arising from reduction in the rho meson mass at finite tempera- 
tures is partially suppressed because of the effect on the width of the 
rho meson induced by Bose enhancement, essentially due to emission 
of pions in a medium of the pion gas in the central rapidity region. 
The sensitivity of the effect on the initial temperature of the hadronic 
phase is also examined. 
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Heavy ion collisions (HIC) provide an unique probe to study the proper- 
ties of hot and dense nuclear matter temporarily produced when two highly 
energetic heavy nuclei collide in the laboratory. Various particles are pro- 
duced directly and via secondary collisions in such high energy HIC which 
serve as a tool to extract information about the environment leading to their 
production; among these the dileptons, as well as the photons, seem to be 
rather promising, as they interact only electromagnetically with the rest of 
the matter and therefore, suffer minimum final state interaction. With this in 
mind several experiments have been performed to measure the dilepton yield 
in HIC m, 0. This could also provide important signals for the novel state 
of matter - the Quark Gluon Plasma (QGP) phase, expected to be produced 
in collisions of two nuclei at ultra-relativistic energies, as per conventional 
wisdom of Quantum Chromodynamics. 

In reality, the contributing factors to this spectra are many and may in 
fact come from the different stages of the evolution of the hot and dense 
hadronic matter temporarily produced in the preliminary stage of the reac- 
tion. In the hadronic stage major contributions may come, say for example, 
from TTTT e~^e~, np vre'^'e", uj 7i^e~^e~, among many others as dis- 
cussed in Ref. In this letter, we demonstrate the effect on low invariant 
mass dilepton yield in heavy ion collisions arising from the increase in width 
of vector mesons (decaying into a pion gas in the central rapidity region) 
due to Bose enhancement (BE) through the induced emission of pions. For 
the sake of illustration and concreteness we study in particular the role of 
the thermal nucleon loop in the rho-mediated production of dileptons, more 
precisely, the dileptons produced from the pion-pion annihilation via Vector 
Meson Dominance (VMD) |Q as shown in Fig. 1. However, the general re- 
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marks regarding the effect of the Bose enhancement would be apphcable to 
any other model. To this end we also provide a 'thumb rule' to estimate the 
extent of the effect of BE at the peak of the dilepton spectrum. 

Actual estimation of the total dilepton spectra requires a thorough under- 
standing of the dynamics of the evolution from the very initial stage of the 
collision. In fact, many calculations, have already been performed, incorpo- 
rating various effects like collision broadening, in-medium effects in the form 
of the effective mass in dense nuclear matter motivated by the novel idea of 
the partial restoration of chiral symmetry, softening of the pion dispersion 
relation [0, H, ||, |^. Still a complete understanding of the data requires fur- 
ther efforts in this direction. With this emerging scenario in the background, 
we calculate the effect of the nucleon loop at finite temperature together 
with the role of the Bose enhancement in the final state taking a particular 
channel into account. It is seen, in this model, that with temperature the p- 
meson mass goes down while on the other hand the p- meson decay width at 
finite temperature shows an interesting behaviour as exposed below. 

The importance of the vrvr e^e~ reaction in the context of high energy 
HIC has already been addressed by many authors |^, especially to explain 
the observed enhancement of the dilepton yield in the low invariant mass 
region as reported by DLS collaboration and most importantly in recent 
times in the SPS 200GeV S+Au reaction. Various attempts have been made 
to account for the data using different models. Li, Ko and Brown explains 
the data by taking the density dependent masses of the vector mesons 
where, however, Bose-enhancement effects on decay width were not taken into 
account. Chanfray observes that medium modified pion dispersion relation 
is responsibile for the excess production of the dileptons in the low mass 
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region ||^, Haglin shows that the other hadronic processes, in particular, 
Tip —>■ ai(1260) 7re+e~, has favourable kinematics to populate masses 
between 2m7r and rrip, although no medium modification for the p meson 
mass or decay width is considered there |T^. In contrast we show that the 
effect of the dropping rho-meson mass is partly offset by the increase in decay 
width of p. 

Theoretical study of tt+tt^ l^l^ dates back to early eighties. Domokos 
and Goldman were among the first to compute the rate for this channel 
using relativistic kinetic theory in which free space pion electromagnetic form 



factor was used and no thermal effects on the p-meson were considered 



Later, Pisarski in the context of chiral symmetry restoring phase transition 
included some finite temperature effects on the p-meson [|^. The VMD 
effects, on the other hand, have been used extensibly at finite temperature 
by Gale and Kapusta, in which, the sensitivity of the dilepton yield on the 
TT — p dynamics in presence of a hot pion gas, has been investigated and they 
concluded that medium corrections are rather modest up to T=150 MeV [|1^]. 
In the present paper, the thermal nucleon loop brings in a modification of the 
rho dominated pion form factor in the medium which is juxtaposed with the 
effect of the thermal distribution of the pions. Unlike the situation in ref. [13 



( concentrated on p7T dynamics ) where the thermal effects through the pion 
loop had little effect, here the nucleon loop at finite temperature is found 
to affect the dilepton yield substantially essentially because the nucleon in 
contrast to the pion suffers a considerable mass reduction. In-medium cross- 
sections are estimated considering both the effect of the thermal distribution 
on the phase space and also its influence on the rho-decay widths in the form 



of the Bose enhancement due to induced emission of pions in a medium of 
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the pion gas. Here a gas of pions are converted into dileptons via rho meson 
by the VMD where necessary temperature dependent effects are included. 
It is observed that the interaction of the rho meson with the thermal bath 
causes an enhancement of the dilepton yield in the low invariant mass region 
because of the lowering of the p- meson mass at finite temperature while 
the Bose enhancement, because of the thermal distribution of pions, causes 
a reduction by increasing the rho decay widths in a hot pion gas. Still, an 
overall excess yield is observed taking both the factors into account together 
with the shift of the peak position towards the low invariant mass region. To 
incorporate the effect of the nucleon thermal loop, we take recourse to the 
usual imaginary time formalism of finite temperature field theory. 
The Lagrangian describing p — N interaction is taken to be 

t"" are the isospin Pauli matrices and N is the nucleon field. The coupling 
constant and the strength of the 'magnetic interaction' is determined either 
from the fitting of the nucleon-nucleon interaction data {gpNN = 2.63, k = 
6.1) or from the VMD of nucleon form factors {Qpnn = 2.71, k = 3.7). 
The form of the polarization tensor for the rho meson is given by 

^M) = E / ^sTrlGik + q, cOn)T,G{k, ^n)f .] (2) 

where /3 = l/kT, T^, T^, are the appropriate vertex factors obtained from the 
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Lagrangian. The sum is carried out over the Matsubara frequencies |T7 , 
ujn = i{'2n + 1)^- This is analogous to what Blaizot discusses in the context 
of photon-self energy for vacuum polarization of the QGP plasma [15]. Here, 
ofcourse we have an additional contribution coming from the tensorial form 
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of the p~N interaction. In the last equation G{k, Un) is the thermal nucleon 
propagator. 



G{k,uOn)= f dre'^-^Gik.T) (3) 

JO 

Gik, r > 0) = A+(A;)7o(l - n^{k))e-'^^ + A_(A;)7on+(A;)e^'=" (4) 



G{k, r < 0) = -A+(A;)7on_(A;)e-^'=" - A-(A;)7o(l - n+(A;))e^'=" (5) 
Here n±{k) = 1/(6^''= ± 1) and A±{k) = ^{^k ± (a ■k + m'jo)), the latter 



is the projection operator as discussed in Ref. |T^ and = ^/W+m*^, 
where the asterisk reminds us that here, instead of the free nucleon mass, 
the medium modified mass is to be used. 

The effective nucleon mass is calculated self-consistently by taking only 
the tadpole diagram arising out of the sigma meson exchange. Here also 



imaginary time formalism is invoked [17] and one subtracts out the divergent 



part as only the shift in the nucleon mass is relevant. 

The covariant rate of the dilepton production for a static thermal system 
is given by the expression 

^ = 1 d%d%fiE,)fiE2)\M\'i27iy5il+ + l_~p,-p2)d%£L (6) 

where d^pi = 2Ei2n)3 /(-^i)? /(-^a) are the pion distribution functions. 
First, we perform the integration over lepton momenta and then the inte- 
gration over the initial momenta are done taking into account the fact that, 
here, the collision is not necessarily colinear. 

The observed dilepton spectra originating from an expanding hadronic 
system (in our case a pionic gas) is obtained by convoluting the static (fixed 
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temperature) emission rate with the expansion dynamics. Relativistic hydro- 
dynamics is a convenient tool to describe the space time evolution of such 



hadronic systems [|I9[. In this work we will use the Bjorken model [^] of 
boost invariant longitudinal expansion, according to which the cooling of the 
system is governed by T(r) = (rj/r)^/^Tj, where Tj is the initial tempera- 
ture at the initial proper time Tj. Within the framework of this model the 
invariant mass distribution of dileptons at mid-rapidity (y=0) is given by, 
dN 3 c dN ^ 1 r , dT 



^ / PTdpT^M%{M)Ko{MT/T) (7) 



dM^dy nRyiak dy ' 2(2? 
where Ra is the radius of the projectile nucleus, c ~ 3.6, = 7r^/30, dN/dy 
is the pion multiplicity, Kq is the zeroth order modified Bessel function of 
second kind. o"j(M) is the in-medium cross section, and also it is to be noted 
that in eq.(6) the matrix element also suffers medium modification because 
of the thermal interaction as represented by the blob in Fig. 1. The in 
medium decay width is also determined in a similar fashion by taking the 
Bose enhancement factor due to the thermal distributions of the pions and, 
in fact, Fp^^TT appears in the denominator of the matrix element M. together 
with the effective nucleon mass. In none of these cases, however, the vertex 
correction is taken into account. 

In Fig. 2, the variation of the effective mass of the nucleon and the p 
meson is presented to show the consistency of our calculation with that of 



others, say for example, Ref. [|T8]. One observes that with temperature (T) 
the p meson mass decreases appreciably which results in a shift of the total 
dilepton yield towards the low invariant mass region. 

The variation of the in-medium decay width (Fp) of the rho meson as 
a function of T which shows an interesting behaviour is depicted in Fig. 
3. Actually the decay width, at finite temperature, is expected to show a 
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reduction as the phase space available is less because of the reduction of the 
p meson mass in the thermal bath, while, on the other hand the presence 
of the pions would cause an enhancement of the p — > tttt decay as evident 
from the Fig. 3. The rapid fall of the p-mass with increasing temperature 
and the consequent reduction in p-width is dramatically offset by the Bose 
enhancement so that upto a temperature of 160 MeV, the width remains in 
the region ~155 to 161 MeV. This particular feature of the p decay width 
has not been remarked upon earlier. The ptttt coupling constant, is estimated 
from the decay width of the rho meson in free space and is given by g^^^/ATr ~ 
2.9. We also note from Fig. 2 that the nucleon effective mass, while showing 
a reduction at finite T, is still large compared to the rho meson mass and 
therefore the nucleon-antinucleon channel would still be closed and have no 
effect on the decay width. 

To bring the effect of Bose enhancement into bold relief, we discuss the 
effect of the nucleon thermal loop on the dilepton spectra with and without 
the Bose enhancement factor included for the p ^ Tin channel. In our ex- 
ample we have taken Ra = 4.6 fm, dN/dy = 225, = 185 MeV, Tp = 130 
MeV, the results are plotted in fig.4, where shift of the peak position towards 
a low invariant mass is observed together with an enhancement of the total 
yield. It is to be noted that inclusion of the Bose enhancement suppresses 
the yield substantially by increasing the decay width as expected from Fig. 
3. Quantitatively, a suppression ~ 30% of the total dilepton yield is observed 
when the effect of BE on decay width is considered compared to the results 
without BE effects everything else remaining the same. However, still an 
excess production of the dileptons at finite temperature is expected once we 
include the effect of the thermal loop. Our findings without the effect of BE 
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on decay width are consistant with what Li, Ko and Brown observes taking 
into account the matter induced modification of the p meson mass. It may be 
noted that upto a temperature(T) ~ 100 — 150 MeV, the BE is dominant as 
evident from Fig. 3 beyond which the mass modification takes over causing 
an over all reduction of the decay width indicating the sensitivity of the total 
dilepton yield on the initial temperature of the hadronic phase. To bring this 
effect into clear focus we also present the spectra for Tj = 165 MeV, keeping 
all the other parameters same as before. The results are displayed in fig. 5. 

The Bose enhancement factor for the rho width due to induced emission 
of pions is clearly given by (l+e"'^"'''/^)^ as each pion carries the energy mp/2. 
Thus the square modulus of the vector meson dominated pion form factor at 

2 

the p-peak (^) receives a reduction factor of (l + e"''^"*''/^)"^ and accordingly 
this would refiect itself on the lepton yield from this channel. However, since 
the basic process has to be integrated over the space-time evolution and the 
concomittant thermal history incorporated in the above reduction factor, an 
effective temperature {T^ff) needs to be inserted. In order to determine T^ff, 
the peak in the invariant mass of the total yield is found and in the given 
model the temperature at which the vector meson mass achieves the value 
TJipeak givcs US an estimate of Te//. Thus the effective reduction factor may 
be estimated to be (1 + e~"*p/^^'^'=//)~^. It is easily checked that this roughly 
agrees with our detailed determination of this value. It should be noted that 
in a case where the vector meson ( say uj ), dominantly decays into three 
pions the suppression due to Bose enhancement will be the sixth power of the 
relevant analogous factor. 

In conclusion, we observe that the thermal nucleon loop and the thermal 
distribution of the pion gas affect the dilepton yield in the low invariant mass 



9 



region quite appreciably and should therefore be included in the dynamical 
calculation to move towards a more complete understanding of the dilepton 
production in high energy heavy ion collision. The main result of the present 
analysis is that the effect of BE on the widths of vector mesons with dom- 
inant pionic decay channels can have substantial influence on the dilepton 
yield depending on the initial temperature. Studies are in progress to inves- 
tigate these effects for other hadronic channels which may contribute to the 
production of leptons. 
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Figure captions 



1. Feynman diagram for dilepton production via Vector Meson Dominance 
from pion-pion annihilation where the blob represents the interaction 
of the p meson with the thermal bath. 

2. Effective masses are depicted as a function of temperature both for 
nucleon and rho meson taking the effect of thermal loops into account. 

3. The variation of the decay width with temperature for three different 
cases are shown. The value of the free decay width is 151 MeV. Solid, 
dashed and dotted curves represent decay width using the free p mass 
with BE, medium modified p width with and without BE respectively. 

4. Effect of the thermal bath on the dilepton yield for = l%bMeV and 
Tp = 130 MeV. Solid, dashed, and dotted represent contribution with 
free (with out medium effect), medium-modified without BE factor and 
medium-modified with BE factor respectively. 

5. Same as Fig.4 for = 165MeV. 
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